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Abstract: The abundance of infrared singularities in gauge theories due to unresolved emission of mass¬ 
less particles (soft and collinear) represents the main difficulty in perturbative calculations. They are typi¬ 
cally regularized in dimensional regularization, and their subtraction is usually achieved independently for 
virtual and real corrections. In this paper, we introduce a new method based on the loop-tree duality (LTD) 
theorem to accomplish the summation over degenerate infrared states directly at the integrand level such that 
the cancellation of the infrared divergences is achieved simultaneously, and apply it to reference examples 
as a proof of concept. Ultraviolet divergences, which are the consequence of the point-like nature of the 
theory, are also reinterpreted physically in this framework. The proposed method opens the intriguing pos¬ 
sibility of carrying out purely four-dimensional implementations of higher-order perturbative calculations at 
next-to-leading order (NLO) and beyond free of soft and final-state collinear subtractions. 
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1 Introduction 

The most important Quantum Field Theories (QFT) for describing elementary particle physics are gauge 
theories. Quantum Yang-Mills theories [1] were introduced about sixty years ago and are now the foundation 
of most of elementary particle theory. One of the main difficulties to extract theoretical predictions from 
perturbative calculations at higher orders is the spread of infinities. In the high-energy regime, there are 
ultraviolet (UV) divergences related to the point-like nature of the theory. Due to the preservation of gauge 
invariance, dimensional regularization (DREG) [2-5], that consists in changing the number of space-time 
dimensions from 4 to d = 4 — 2e, is the customary regularization method to overcome this problem. The UV 
singularities in four dimensions manifest as poles in the dimensional parameter e, and are easily removed 
by renormalization. The real challenge is, however, in the infrared (IR) regime. Unresolved radiation of 
massless particles (soft and collinear) is physically degenerate with the absence of emission, while from the 
theory point of view they are described as different final states. This leads to further infinities that are also 
regularized in DREG. 

The standard approach [6-14] to remove the IR singularities consists of adding to the real radiation 
contribution suitable subtraction terms that mimic the IR behavior of the emission scattering amplitudes and 
render them finite in the IR limit. The same subtraction quantities, integrated over the phase-space of the 
extra radiation, are subtracted back from the virtual corrections. The sum of both contributions leads to finite 
theoretical predictions for physical observables, if all the degenerate states are considered in the sum. In this 
paper, we describe an alternative approach which is based on the loop-tree duality (LTD) theorem [15- 
22]. We propose a new method that combines simultaneously virtual and real corrections, following in 
fact the original aim of the LTD method. We present here the first practical realization of LTD that we 
illustrate with some reference examples. This fact has also a strong implication: the possibility of carrying 
out purely four-dimensional implementations free of soft and final-state collinear subtractions. Unlike the 
method proposed in Ref. [23], LTD does not require any smearing in the total energy to combine real and 
virtual corrections. Likewise, massless internal propagators are kept massless in LTD, as opposed to the four 
dimensional implementations presented in Refs. [24-26]. 




2 Review of the loop-tree duality 


The LTD theorem [15] establishes a direct dual relation between one-loop 7V-leg scalar integrals or scattering 
amplitudes in any relativistic, local and unitary QFT with corresponding tree-level objects integrated in a 
modified phase-space that resemble real corrections. Specifically, the dual representation of a A^-leg scalar 
one-loop integral consists of the sum of N dual integrals: 
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( 2 . 1 ) 


where GD{qi',qj) = {q] ~ ~ dual propagators, with i,j £ ai = {1,2,.. .N}. The 

momenta of the internal lines are denoted = {qi^, q^), where o is the energy (time component) and 

are the spacial components. They are defined as g^ = i+ki with £ the loop momentum and ki = pi + .. .+Pi. 
The four-momenta of the external legs are pi, which are taken as outgoing, with kN = 0 by momentum 
conservation. The loop integration measure reads 
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The delta function 5 (qi) = 27r i 0(gi o) ^{ql — 'm'i) sets the internal lines on-shell and is the result of selecting 
the poles of the Feynman propagators with positive energy g^ o and negative imaginary part by applying the 
Cauchy residue theorem along a contour that is closed in the lower half of the complex plane. LTD is 
equivalent to perform the loop integration in the d — 1 dimensional space defined by the forward on-shell 
hyperboloids, Gi?(gi)“^ = (gf—mf-|-i0) = Oandgio > 0, or forward light-cones for massless propagators. 

The key feature of the dual representation in Eq. (2.1) is that the usual Feynman propagators have been 
replaced by dual propagators whose -fiO prescription depends now on the sign of p • kji, where rj is a future- 
like vector, > 0, with positive definite energy po > 0, and kji = qj — qi. The idea of having different 
-fiO prescriptions for different propagators is at first sight astonishing, but indeed it is a necessary condition 
for the consistency of the method. As shown in Ref. [18], the integrand in Eq. (2.1) becomes singular at 
the intersection of forward on-shell hyperboloids (light-cones for massless propagators), and forward with 
backward (g^ o < 0) on-shell hyperboloids. Those singularities lead to integrable threshold singularities 
or non-integrable soft and collinear divergences. In the forward-forward case, there is a cancellation of 
singularities among different dual contributions, and the change of sign of the dual prescription plays a 
central role in that cancellation. In the forward-backward case, the singularities remain constrained to a 
compact region of the loop three-momentum and admit a nice physical interpretation in terms of causality, 
which is indeed consistent with the Cutkosky rule. These singularities occur only in one direction, the 
direction of the chosen internal momentum flow, when the on-shell virtual particle is emitted and interacts 
with another outgoing virtual on-shell particle (gj o < 0) after loosing energy by radiating external particles. 


3 Dual cancellations 

To validate the LTD method with an explicit example, we focus our attention to the simplest case of the 
scalar three-point function, with pa ^ pi + P2 as external momenta and pg = si2 > 0; pi, p2 and the 
internal momenta are taken massless. The loop integral is defined by the internal momenta qi = £ + pi, 
q 2 = £ + Pi 2 and gg = £, according to Eig. 1. We proceed step by step, by integrating first in DREG and 
then by repeating the calculation in a suitable representation where the limit e —0 can safely be taken at the 
integrand level. The dual representation of the scalar three-point function consists of three contributions: 

h = - fHq^)T\GD{q^■,qJ), z,j = l,2,3. (3.1) 
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Figure 1. Kinematic configuration of the three-point function for the process pa ^ pi + P2, its interference with the 
Bom amplitude (left), and interference of real radiation contrihutions with ps —^ p'l + P2 + p'r (right)- 


Assuming pi oriented along the positive z-axis, p 2 in the opposite direction, the on-shell loop momenta are 
parametrized as 

^ (1, 2\/tt,(l - V,) e,,i, 1 - 2vi) , (3.2) 

with 0 S [0, 00 ), Vi G [0,1], and a unit vector in the transverse space. Taking into account these 
parametrizations the scalar products are given by 2 qi -pxjsyi = Ci,o Vi, and 2 qi ■ P 2 /S 12 = Ct.o (1 ~ Vi), and 
the dual integrals are written as 
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h 
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where we defined the d-dimensional integration measures as 

512 Si,0 I 

d[vi] = (t;i(l - Vi))~^ dvi , (3.4) 


and the dual vector was taken as = (1, 0). Figure 2 (left) shows a graphical representation of the cor¬ 
responding light-cones and their intersections where threshold and IR singularities appear. In Eq. (3.3), 
only the dual prescription of the dual integral I 2 has been kept explicitly, as it regulates a threshold singu¬ 
larity. The dual prescription of the other propagators is required for the consistent partial cancellation of 
singularities among dual integrals [18], but can be removed for this practical calculation. 

The explicit computation in DREG of the total integrals in Eq. (3.3) gives 

Ii = 0 , I 2 = cr ^ I I 3 = cr ^ ) (3-5) 


with 

.. _ 1 r(l + e)r^(l-e) 

cos(7re) ’ (47r)2-'^r(l — 2e) 


(3.6) 


where cr is the volume factor of the phase-space integrals, and cr is the volume factor of the loop integrals. 
The dual integral Ii vanishes because it is scaleless, and following the usual convention in DREG it is 
analytically continued by matching its UV and IR poles, i.e. euv = cir- Summing up the results in Eq. (3.5) 
we obviously obtain the well-known result 


« o 2e 
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Figure 2 . Light-cones of the massless three-point function in the loop coordinates I = ^812 !“2 {^0, ^±, ^z), with 
?-L = \/Q + (left). Each forward light-cone correspond to the integration region of one dual integral U, and 
the intersection with the other light-cones generates the singularities of the integrand. Location of threshold and IR 
singularities in the loop three-momentum space (right). 


At this point, we start the discussion by isolating the singular regions of the loop integrand, which are 
bounded to a compact region of the loop three-momentum, as it can be seen from Fig. 2 (right). The dual 
integral Ii develops a soft singularity at ^10 = 0 and a collinear singularity at ui =0 with ^10 S [0,1], 
namely G [—1,0]; the collinear singularity at 6 > 0 cancels with the other dual contributions. Likewise, 
the dual integral I 2 is collinear at U 2 = 1 with ^ 2,0 G [0jl]> or fz G [~li0]> ^rid undergo cancelling 
collinear divergences at 6 < -1- The integrand of I 2 gets in addition a threshold singularity at ^ 2,0 = 
which prevents us from having a well defined integral if we restrict the integration region by ^ 2,0 = 1- 
Consequently, we introduce an energy-cut w, with w close to zero, and define the following soft and collinear 
contributions to the dual integrals 

= /i(a.o < w) , 

=/i(u;<ei.o<l;t^i<l/2), 

= /2(6.o < 1 + > 1/2) . (3.8) 

The soft integral ^ contains the soft singularity of the dual integral Ii, and its collinear singularities close 
to the soft region, ^10 < w. Once the soft part is removed, the remaining collinear singularities of Ii 

(c) 

are contained in the collinear integral I) . It is necessary to extend the integration limits of the collinear 

(c) 

integral beyond the strict collinear singular region of I 2 in order to treat properly its threshold singularity 
at ^ 2.0 = 1- Integrating these contributions in DREG, we obtain the following results: 
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with 2 T 1 (2e, 2e, 1 -|- 2e; z) = 1 -|- 4e^ Li 2 ( 2 :) + 0{e^). As expected, the soft integral contains double 

fc) (c) 

poles in e, while the collinear integrals J) and develop single poles only. Each integral depends on the 
cut w, but the IR poles of the sum are independent of w and they agree with the total divergences of the full 
integral in Eq. (3.7). 
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At this stage, we have shown for the first time with a concrete example that the IR singularities in the 
LTD representation are constrained to a compact region of the loop three-momentum, as it was anticipated 
before based on the general demonstration from Ref. [18], and have recalculated the IR poles only by taking 
into account this compact region. The definition of the soft and collinear integrals is, however, somehow 
arbitrary; strictly, the IR singularities of the loop integrand are constrained to the region defined by the limit 
w 0 and more precisely to the segment defined by € [—1,0]. Any compact region of the loop three- 
momentum enclosing € [—1,0] will lead to the same e poles, although to different finite contributions. 
We will later redefine the IR regions of the loop integrand in a more convenient way to combine them with 
the real corrections. 

Outside the soft and collinear regions, the sum of the dual integrals is finite because the integrand 
singularities appear at the intersection of forward light-cones. We define first a forward region with Vi < 1/2: 
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The dual integrals Ii and I 3 are divergent at ui = U 3 = 0 in the forward region, but their sum is finite in this 
collinear limit. However, they diverge in the UV in all directions in spite of the absence of UV divergences 
in the original integral. This is due to the fact that dual propagators are linear in the loop momentum. The 
UV finiteness of the original integral is recovered by adding together all the dual integrals [15]. Therefore, to 
cancel this UV behaviour it is necessary to add the contribution of l 2 - The integral over the sum of the three 
dual integrands in Eq. (3.10) can thus be performed with e = 0. In this expression, we have also identified 
all the angular and energy variables, Vi = v, ~ 1 = '^ 2,0 = ^ 3,0 = ■Co. to obtain common integration 
limits. Notice that it is necessary to keep the dual iQ prescription of I 2 because its threshold singularity at 
^ 2,0 = 1 appears within the integration region; it does not affect the UV cancellation. The logarithmic term 
of the integrand in Eq. (3.10) is the result of the mismatch in the e-expansion of the integration measure of Ii 
and lo- The same result is obtained without extra logarithmic contributions by expressing (^ 1 , 0 , ui) in terms 
of (^ 3 , 0 ; ^^ 3 ) at the price, however, of introducing cumbersome integration limits. We obtain in both cases 


/(f) =cr — 
S12 



in log(2) 


+ 0 {e). 


(3.11) 


The other finite contribution in the backward region (vi > 1/2) is 


/(f’) = fl(fl,0 > W -jVi > -) + 12(^2 ,0 > 1 + W ; W2 > 2) “f -^3(^3 > 2^ 


= Cr 


S12 


2 Li 2 ( -— ) - log(2) log(u;) 
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(3.12) 


Now, the collinear cancellation occurs at ui = 1 = t ;2 and the UV behaviour of Ii + I 2 is cancelled by 
I 3 . As expected, the sum of the soft, collinear, forward and backward integrals in Eq. (3.9), Eq. (3.11) and 
Eq. (3.12) leads to the full result given in Eq. (3.7), up to 0{e). 


4 Real-Virtual mapping and cancellation of infrared singularities 

In the previous section we have analyzed the divergent structure of a scalar one-loop integral, and have 
isolated its IR behaviour to a compact region of the loop three-momentum. We shall illustrate now how 
to combine virtual corrections with real radiation in order to obtain physical cross-sections. To show the 
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advantages of the method it is enough to work out a simplified calculation where the one-loop scattering 
amplitude is just given by the scalar three-point function 


= -ig^ si2L^^\pi,P2,-P3) , 


(4.1) 


where g is an arbitrary coupling. Explicit computations in a real physical model will be presented elsewhere 
[29, 30]. The corresponding scalar tree-level amplitude is {pi,p 2 ',P 3 )) = ig, and we shall take into 
account their interference. Re integrated over the two-body phase-space. We consider also 

scalar scattering amplitudes with emission of an extra particle \M\^\p[,P 2 ,Pr]P 3 )) = —ig^y/^ls'^r > 
in particular the interference (see Fig. 1 (right)) 


Re{MZ’\MfJ)=g 


(4.2) 


which is integrated over the corresponding three-body phase-space ^ . To combine virtual and real correc¬ 
tions directly, we first decompose the real-emission phase-space by using the identity 1 = 0 {y 2 ^ — y[^) + 
9{y'ir — 2 / 2 t-)> where = s'^/si 2 . It is unnecessary to single out a soft region. In general, segmenting the 
real-emission phase-space into different regions featuring at most one soft and/or collinear singularity is a 
common practice in many subtraction schemes, as for instance in FKS [ 8 ]. The main novelty of the FTD ap¬ 
proach is that that segmentation will be translated into the virtual contributions such that a local cancellation 
of singularities at the integrand level is achieved. 

Thus, we define the following dual contributions to the total cross-section 

= ^/ d^i^3 2Re{Mt^\MfJ)0iy'^,-yl), 

= ^ / 2Re |7Wf)) e{y'^, - y^) , i, j = 1, 2 , (4.3) 

where \Ai^P) = —ig^ Si 2 h is the dual component of the one-loop scattering amplitude. The third dual 
component given by I 3 in Eq. (3.3) does not contribute to the IR structure of and can be integrated 

independently of the real corrections. 

The overall picture presented in Eq. (4.3) must be completed with a mapping of the virtual and the real 
kinematics. For a fixed value of the external momenta pi and p 2 entering the loop scattering amplitude, and 
the loop three-momentum £, we define the following mapping, with qi = i + pi on-shell. 


pT = Ql : P'l = -03 + OilP2 = Pi - 01 + «1 P 2 : 

P 2 = (1 - ai)P 2 ) «1 = 77 ^— ) (4.4) 

2^3 • P2 

where qs = qi—pi is off-shell. Momentum conservation is automatically fulfilled as p[ +P 2 +p'^ = pi +P 2 , 
and all the primed final-state momenta are on-shell. The non-collinear parton P 2 plays here the role of the 
spectator parton in the dipole formalism [ 6 ] . The mapping in Eq. (4.4) is also valid for an arbitrary number of 
external particles, provided that p' = pi for partons not entering in the collinear subprocess. The two-body 
invariants are given by 


Vlr = 


Vl 6,0 


1 - (1 - 


/ (1 - fi) (1 - 6 . 0 ) 6.0 

2/2r = 


1 - (1 - 


2/12 = 1 - 6.0 


(4.5) 


These expressions are obtained from the momentum parametrizations given in Eq. (3.2) and Eq. (4.4). The 
loop momentum qi andpi become collinear at v\ 0, with ai —0 and 0 in this limit. According to 

^For identical final-state particles all the other possible interferences should also be considered. Besides, the squared amplitudes 
are necessary, which have to be mapped to self-energy insertions of external legs. For the current illustrative purposes of the 
paper, however, it is enough to consider this interference. 
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Eq. (4.3), we use this mapping in the regions of the two-body and three-body phase-space where y[j. < y^r- 


= ^v,i + ^R,i = 25 ^ J c^Ki.o] d[vi] 0(1 - 2i;i) 0 

1 — '^ 1,0 


1 - 2 ui 


Vi 


-Ci.t 




1 - (1 - Ui)^i,| 


- 1 


(4.6) 


where = y^/(2si2) / d$i _>.2 would be the lowest order cross-section. In Eq. (4.6), the three-body 
phase-space d$i _>.3 was written in terms of the loop variables ^10 and vi by using Eq. (4.5). The integrand in 
Eq. (4.6) has the form f{vi, ^ 1 , 0 ), but f{vi = 0) = 0 = /(^i,o = 0) with f{vi, ^ 1 , 0 ) = 0{e), 

and thus this integral is of order e, i.e. = 0(e). Similarly, we dehne a complementary mapping to cover 
the remaining three-body phase-space, which includes the other collinear limit with ^ 0. Considering 
52 on-shell, the corresponding mapping is 


/a a 

P 2 = 92 : 


'U Ml M M M I M 

Pr = -9l + a2Pl = P 2 - 02 + «2Pi , 


Pi = (1 - “ 2 ) Pi 

and the associated invariants are 
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1 — U 2 ^ 2,0 
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— '^ 2,0 ) (1 — V 2 ) 


-1 


(1 — V 2 ^ 2 , 0 )^ 


1 ~ '^ 2,1 


■iQ 


- ITlSil - ^ 2 , 0 ) 


and we obtain directly with e = 0 


6 


0 {e), 


(4.9) 


(4.10) 


where a = / (dir)^. The condition < y[j. excludes the threshold singularity of I 2 from the integration 

region with the exception of the single point at V 2 = 1, C 2.0 = 1- This makes unnecessary the introduction 
of an energy-cut w as in Sec. 3. 

Einally, analogously to the forward and backward dual integrals defined in Sec. 3, we collect the rem¬ 
nants of the three dual integrals obtained after excluding the respective regions of the loop three-momentum 
already included in Eq. (4.6) and Eq. (4.10). Explicitly, we dehne 


ff!;> = »<»> 


29 " 


- J d[^i,o]d[vi] ^1 - 0(1 - 2ui) 0 ^ 

— J ^[^ 2 . 0 ] ^^ 2 ] 0 ^^2,0- ^ 

+ J f^[6,o] ((ba] 


1 - 2 ui 
1 — Ul 

( 1 -U 2 )”^ 


“ ^1.0 ) ) Ci,o ('*^1(1 “ '*^1)) 
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-1 


1 ~ '? 2,0 


iO 


+ 171(5(1 - 6,0)^ 


(4.11) 


Although the individual contributions in Eq. (4. 1 1) are singular, this virtual cross-section is UV and IR hnite, 
and can be calculated with e = 0. We apply the change of variables used to obtain Eq. (3.11) and Eq. (3.12) 
to avoid logarithmic contributions from the e-expansion of the integration measure. The result is 




(4.12) 
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The sum of all the contributions, Eq. (4.6), Eq. (4.10) and Eq. (4.12), gives a total cross-section of 0{e), in 
agreement with the result that would be obtained from the standard calculation in DREG. 

To conclude this section, we want to emphasize that DREG was only used for a consistency check of 
our results. The advantage of the method presented here is that NLO cross-sections can be computed in 
four dimensions, i.e. with e = 0, because soft and collinear divergences of virtual and real corrections are 
matched locally at the integrand level. Collinear factorization guarantees that this matching is fulfilled in 
QCD [18, 27]. Threshold singularities should be treated apart in numerical calculations, e.g. by contour 
deformation [23, 28]. A first numerical implementation of multi-leg scalar and tensor one-loop integrals in 
LTD has been presented in Ref. [19-21]. 

5 Ultraviolet Renormalization 

LTD also offers an appealing physical interpretation for the renormalization of UV divergences at the inte¬ 
grand level. To illustrate the cancellation of these divergences, we consider the simplest purely UV scalar 
two-point function, with massless internal lines. A detailed calculation of this scalar two-point function in 
the LTD framework has been presented in Ref. [15]. Eollowing the suggestion of Ref. [28], we construct its 
UV counter-term from 

rent _ 

— 

where guv = £ + kuy, with fcuv located somewhere in the loop momentum space. This counter-term has 
the same local UV behaviour as the integrand of the scalar two-point function in all the loop directions. We 
shall now find the dual representation of I{jy in three dimensions. 


I 


^ (^uv 


(5.1) 


^0 



Figure 3. On-shell hyperboloids of the ultraviolet counter-term. 


The dual representation of Eq. (5.1) requires to deal with double poles in the loop energy [17], with the 
UV propagator being massive. The calculation of the corresponding residue leads to the dual representation 


rent _ 

''UV — 


I 


s (guv) 



(5.2) 


where g^y g = Quv + Muv ~ detailed calculation of the double pole residue has been presented 

in Ref. [17]; again we used = (1, 0) as dual vector. The arbitrary scale /ruv admits a direct physical inter¬ 
pretation as renormalization scale because the UV counter-term only contributes for loop energies larger than 
—kuv.o + Muv in the LTD representation (see Eig. 3), although with unconstrained loop three-momentum. 
As renormalization scale, p,uv has to be chosen of the order of the physical hard scale. In addition, a clear 
criteria to select /iuv and kjjv arises if we pretend that the on-shell hyperboloids of the UV propagator 












in Eq. (5.1) do not intersect with any of the forward on-shell hyperboloids of the original integral. Since the 
distance between the UV forward and backward on-shell hyperboloids is 2/iuv> the minimal choice fulfill¬ 
ing the aforementioned conditions would be 1 /2 of the hard scale, which accidentally agrees with one of the 
standard limits used to estimate the theoretical uncertainty. Of course, this is an ad hoc argument that is not 
supported by any physical evidence, but it represents a good interpretation of that choice. 

In consequence, subtracting the UV counter-term from the scalar two-point function we obtain 


lW(p,-p)-^uv = 


(47r)2 


- log - 




- iO -f 2 


0{e), 


(5.3) 


and this calculation gives the same result if each of the individual contributions is integrated in DREG or if 
the sum of the integrands is taken at e = 0. 

We conclude this section by noting that contrary to the four-dimensional UV renormalization introduced 
in Ref. [24], where the limit fi ^ 0 is taken with the purpose of removing the polynomial dependence in fi 
and the remainder is evaluated at the renormalization scale fi = /tr, the counter-term defined by Eq. (5.1) 
keeps the complete dependence in puv- Eventual polynomial dependences in puv can be cancelled by 
adding to the counter-term UV subleading contributions such that only the UV pole is subtracted [28]. 


6 Conclusions 

We have applied the LTD formalism to physically reinterpret perturbative calculations in QET. With some 
reference examples, we have reanalyzed the infrared and causal structure of NLO computations. The sum 
over degenerate IR states is accomplished at integrand level by mapping the phase-space of real radiation 
with that of the virtual contributions in the region of the loop three-momentum which is of the size of the 
physical scales. This represents a new paradigm in perturbative calculations as virtual and real corrections are 
not regularized independently, but they act directly as the IR subtraction counter-term of each other. Outside 
this region, at large loop three-momentum, virtual corrections require only the subtraction of UV divergences 
at integrand level to render them finite and renormalization acquires a clear physical interpretation. 

The application of LTD to loop scattering amplitudes in gauge theories was already discussed in Ref. [15], 
and the momentum mappings introduced in this paper can be extended straightforwardly to massless n-body 
scattering processes, and even generalized to the massive case. Each mapping is suitable for a different 
collinear configuration, and has to be applied to the sum of diagrams belonging to gauge invariant subsets 
featuring the same collinear limit. 

The results presented in this paper represent a promising first effort towards a four-dimensional imple¬ 
mentation for the computation of physical cross-sections at NLO without introducing soft and final-state 
collinear subtractions. They also open the attractive possibility of extending this approach to next-to-next-to 
leading order (NNLO) and beyond. Eurther details will be given in a forthcoming publication, including the 
application to gauge theories [29]. 
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